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Abstract 

An improved measurement of the process —* p^tr is presented. The data sample of 211 fb _1 
comprises 232 million T(4S) — > BB decays collected with the BABAR detector at the PEP- 
II asymmetric B Factory at SLAC. The yield and CP asymmetry are calculated using an ex- 
tended maximum likelihood fitting method. The branching fraction and asymmetry are found 
to be B(B ± -> p ± 7r°) = [10.0 ± 1.4 (Stat.) ± 0.9 (Syst.)} x lO" 6 and A C p(B ± -> ^vr ) = 
—0.01 ± 0.13 (Stat.) ± 0.02 (Syst.), superseding previous measurements. The statistical signifi- 
cance of the signal is calculated to be 8.7a. 
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1 Introduction 



Branching fraction and CP asymmetry measurements of charmless B meson decays provide valuable 
constraints for the determination of the unitarity triangle constructed from elements of the Cabibbo- 
Kobayashi-Maskawa quark-mixing matrix ^H]. I n particular, the angle a = arg [— V^VJ^/V hY^A 
of the unitarity triangle can be extracted from decays of the B meson to p^ir^ 1 final states [3]. 
However, the extraction is complicated by the interference of decay amplitudes with differing weak 
and strong phases. One strategy to overcome this problem is to perform an SU(2) analysis that uses 
all pn final states [3] . Assuming isospin symmetry, the angle a can be determined free of hadronic 
uncertainties from a pentagon relation formed in the complex plane by the five decay amplitudes 
B° p+TT~, B° -> p~Tr + , B° -> p vr°, B + p + ir° and B + -> p°ir + . These amplitudes can 
be determined from measurements of the corresponding decay rates and CP asymmetries. While 
all these modes have been measured, the current experimental uncertainties need to be reduced 
substantially for a determination of a. Here we present an update to a previous measurement jS] 
of the B^ — ► p^ir branching fraction and CP asymmetry 

N(B~ -» p--K°) - N(B + -» p+vr°) 
CP ~ N{B~ p-ir°) + N(B+ ^ p+7T°) ' 



2 Data Set and Candidate Selection 

The data used in this analysis were collected with the BABAR detector 6 at the PEP-II asymmetric- 
energy e + e~ storage ring at SLAC. Charged-particle trajectories are measured by a five-layer 
double-sided silicon vertex tracker and a 40-layer drift chamber located within a 1.5-T solenoidal 
magnetic field. Charged hadrons are identified by combining energy-loss information from tracking 
with the measurements from a ring-imaging Cherenkov detector. Photons are detected by a CsI(Tl) 
crystal electromagnetic calorimeter with an energy resolution oia E /E = O.O23(S/GeV)- 1 /4 0O .oi4. 
The magnetic flux return is instrumented for muon and K® identification. 

The data sample includes 232 ± 3 million BB pairs collected at the T(4f>) resonance, corre- 
sponding to an integrated luminosity of 211 fb . It is assumed that neutral and charged B meson 
pairs are produced in equal numbers [7]. In addition, 22 fb -1 of data collected 40 MeV below the 
T(4S) resonance mass are used for background studies. 

We perform full detector Monte Carlo (MC) simulations equivalent to 460 fb -1 of generic BB 
decays and 140 fb -1 of continuum quark-antiquark production events. In addition, we simulate 
over 50 exclusive charmless B decay modes, including 1.4 million signal B^ — ► p^ir decays. 

B meson candidates are reconstructed from one charged track and two neutral pions, with the 
following requirements: 

Track quality. The charged track used to form the B^ — > p^ir candidate is required to have 
at least 12 hits in the drift chamber, to have a transverse momentum greater than 0.1 GeV/c, and 
to be consistent with originating from a S-meson decay. Its signal in the tracking and Cherenkov 
detectors is required to be consistent with that of a pion. We remove tracks that pass electron 
selection criteria based on dE/dx and calorimeter information. 

7T° quality. Neutral pion candidates are formed from two photons, each with a minimum energy 
of 0.03 GeV and a lateral moment [S] of their shower energy deposition greater than zero and less 
than 0.6. The angular acceptance of photons is restricted to exclude parts of the calorimeter where 
showers are not fully contained. We require the photon clusters forming the 7r° to be separated in 
space, with a ir° energy of at least 0.2 GeV and an invariant mass between 0.10 and 0.16 GeV/c 2 . 
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Kinematic requirements. Two kinematic variables, AE = E* B — yfsjl and the beam energy 



substituted mass of the S-meson toes = y (s/2 + p • pb) 2 /Eq — p B , are used for the final selection 

of events. Here E* B is the energy of the B meson candidate in the center-of-mass frame, Eq and \fs 
are the total energies of the e + e~ system in the laboratory and center-of-mass frames, respectively; 
Po and pb are the three-momenta of the e + e~ system and the B candidate in the laboratory frame, 
respectively. For correctly reconstructed p^n candidates AE peaks at zero, while final states with 
a charged kaon, such as B^ — > K^ir , shift AE by approximately 80 MeV on average. Events 
are selected with 5.20 < toes < 5.29 GeV/c 2 and \AE\ < 0.20 GeV. The AE limits help remove 
background from two- and four-body B decays at a small cost to signal efficiency. 

Continuum suppression. Continuum quark-antiquark production is the dominant back- 
ground. To suppress it, we select only those events where the angle #g ph in the center-of-mass 
frame between the direction of the B meson candidate and the sphericity axis of the rest of the 
event satisfies |cos#g ph | < 0.9. In addition, we construct a non- linear discriminant, implemented 
as an artificial neural network, that uses three input parameters: the zeroth- and second-order Leg- 
endre event shape polynomials Lq, Li calculated from the momenta and polar angles of all charged 
particle and photon candidates not associated with the B meson candidate, and the output of a 
multivariate, non-linear B meson candidate tagging algorithm jjjj. The output ANN of the artifi- 
cial neural network is peaked at 0.5 for continuum-like events and at 1.0 for B decays. We require 
ANN > 0.63 for our event selection. 

p mass window. To further improve the signal-to-background ratio we restrict the invariant 
mass of the p candidate to 0.55 < m W7T < 0.95 GeV/c 2 . 

Multiple candidates. Neutral pion combinatorics can lead to more than one I?-meson can- 
didate per event. We choose the best candidate based on a \ 2 formed from the measured masses 
of the two 7T° candidates within the event compared to the known ir° mass |10| . In the case of 
multiple charged pion candidates the choice is random. 

Efficiency. The total B^ — * / 9 ± 7r° selection efficiency is 15.4 ± 0.1%. In MC studies, the 
signal candidate is correctly reconstructed 54.9% of the time. The remaining candidates come from 
self-cross-feed (SCF, 37.5%) and mistag events (7.6%). We define SCF events as those where one 
or more elements of the l?-candidate reconstruction are incorrect except for its charge. They stem 
primarily from swapping the low energy tt° from the resonance with another from the rest of the 
event. Signal events reconstructed with the wrong charge are classified as mistag events. Both 
SCF and mistag events emulate signal events, however the resolution in toes an d AE tends to be 
worse. 

3 Background Contributions 

MC events are used to study backgrounds from other 5-meson decays. The dominant contribution 
comes from b — > c transitions; the next most important is from charmless £>-meson decays. The 
latter tend to be more problematic as the branching fractions are often poorly known, and because 
they may peak at the same invariant mass as the signal B^ — ► p^ir events. Seventeen individual 
charmless modes show a significant contribution once the event selection has been applied (Table^). 
These modes are added into the fit fixed at the yield and asymmetry determined by the simulation, 
given an assumed branching fraction. Wherever branching fractions are not available, we use half 
the upper limit. If no charge asymmetry measurement is available, we assume zero asymmetry. 
p* resonances. Although all other states which decay like the p to tttt - subsequently referred 
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Backgrounds to 


p^n from S-related Sources 




Mode 


Efficiency 

(%) 


Assumed B 
(xHT 6 ) 


Assumed ,A.cp 


Expected 
yield 


Generic b — > c 


1.69 x 10" 4 


- 


0.00 


392.3 ± 19.8 




P ± P T 


2.12 


30.0 ±6.0 HI] 


0.00 ± 0.00 


147.7 ± 29.5 


_B° — > 




2.39 


24 ± 2.5 dU 


0.00 ± 0.20 


132.9 ± 13.8 


B±- 




4.12 


5.5 ±0.6 HI] 


-0.02 ± 0.07[II] 


52.5 ±5.7 


S ^ 


a?7r° 


1.16 


17.5 ± 17.5 


0.00 ± 0.20 


46.9 ± 46.9 


B ± - 


p±p° 


0.57 


26.4 ±6.4 [IT] 


-0.09 ±0.16 [II] 


34.6 ± 8.4 


B ± - 


"K^~K S [ — > 7r°7r°) 


1.49 


3.74 ± 0.20 [II] 


-0.02 ± 0.03QT] 


12.9 ±0.7 


B ± - 




0.41 


12.1 ±0.8 mi 


0.04 ± 0.04 mi 


11.5 ±0.8 


B ^ 


7r°7T 


2.62 


1.51 ±0.28 [II] 


0.00 ± 0.00 


9.2 ± 1.7 


B° — > 


t/tt 


1.15 


1.85 ± 1.85 jll] 


0.00 ± 0.20 


4.9 ± 4.9 


S±- 




0.17 


9.76 ± 1.22 [II] 


0.00 ± 0.20 


3.8 ±0.5 




TT°K*° 


0.90 


1.7 ±0.8 mi 


0.00 ± 0.20 


3.5 ± 1.7 


s±- 




0.61 


2.3 ±0.8 [12] 


0.04 ± 0.29 


3.4 ± 1.0 




il-ipT 


0.14 


9.9 ± 1.6 mi 


0.17±0.16HU 


3.2 ±0.5 


B±- 




0.02 


40.3 ±2.6 [HI 


-0.01 ± 0.03HH 


2.1 ±0.1 


S±- 


P ± 7 


0.65 


0.9 ±0.9 mi 


0.00 ± 0.20 


1.4 ± 1.4 




K^pT 


0.05 


12.0 ± 12.0 [II] 


0.00 ± 0.20 


1.4 ± 1.4 


s±- 


p ± (1450)vr° 






0.00 ± 0.20 


8±8 


Total 








872.2 ±62.1 



Table 1: The individual S-background modes considered in the fit. The expected number of events 
after all cuts are listed for each mode. For modes which do not have well measured branching 
fractions, half the upper limit is used. If no Acp measurement is available, we assume zero 
asymmetry with a 20% uncertainty. 
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to as p* - lie outside our p(770) mass cut, a contribution to our signal cannot be a priori ruled 
out. The only non-strange vector resonances which can decay to two pions are the p(1450) and 
the p(1700). To account for the possible presence of these modes, a fit to the — > / o* ± 7r° yield is 
performed in a sideband of the invariant mass using the three variables mES> A.E, and ANN. The 
mass window is chosen to be as far as possible from the p(770) mass, centered near the pole of the 
,9(1700) at 1.5 < m niT < 2.0 GeV/c 2 . The fitted yield for the B ± — > / 9* ± 7r° decay is then extrapolated 
into the nominal region. Although the choice of mass range is motivated by the /a(1700), any yield 
seen is attributed entirely to the /o(1450), which is the closer of the two resonances to the signal. 
From the — > /? ± (1450)7r° MC, the ratio of candidates in the sideband to candidates in the signal 
mass region is approximately 12.6:1. The fit in the sideband yields 101 ± 32 events, resulting in 
an estimate of the p* background of 8 events. We assign a conservative systematic uncertainty 
of 100% to this number. The p* then enters into the nominal fit with PDFs constructed from 
B ± -► ^(1450)71-° MC simulation. 

Non-resonant decays to 7r ± 7r°7r°. The non-resonant B^ — > tt^tt^tt branching fraction has, 
to date, not been measured. To estimate the significance of its contribution we select a region of the 
Dalitz plot — defined by the triangle (m^.m^o) = (6,6), (6, 15), (11, 11) GeV 2 /c 4 — that is 

far from the signal as well as p(1450) and higher resonances and which has low levels of continuum 
background. A likelihood fit in this region yields —5.1 ± 7.6 non-resonant events in a data sample 
of 1100 events. This is consistent with zero. The non-resonant contribution is therefore deemed 
negligible. 



4 The Maximum Likelihood Fit 

An unbinned maximum likelihood fit to the variables ?ties> A.E, m^, and ANN is used to extract 
the total number of signal B^ — > p^tr and continuum background events and their respective 
charge asymmetries. The likelihood for the selected sample is given by the product of the probability 
density functions (PDF) for each individual candidate, multiplied by the Poisson factor: 

1=1 

where N and N' are the number of observed and expected events, respectively. The PDF Vi for a 
given event % is the sum of the signal and background terms: 

Vi = iV Sig xi[(l-Q^ Si s)/ Sig Pf s 

+ (l-Q^ Si S)/sCF^ g Fii 

+ (l + QiA s ^)f Mis V^.} 
+ ^Nf k ^x 1 -(l-Q t A^)vf^, 

3 

where Qi is the charge of the pion in the event, iV Si s(A^ Bkg ) and A Si z(Af kg ) are the yield and 
asymmetry for signal and background component j, respectively. The fractions of true signal 
(/sig)> SCF signal (/scf), and wrong-charge mistag events (/Mis) are fixed to the numbers obtained 
from MC simulations (Section |5J). The j individual background terms comprise continuum, b — > c 
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decays, p* , and 17 exclusive charmless B decay modes. The PDF for each component, in turn, 
is the product of the PDFs for each of the fit input variables, V = Vm ES AE'PANN'Pm, 7T7T - Due to 
correlations between AE and niES, the "P mE s,A£: for signal and all background from B decays are 
described by two-dimensional non-parametric PDFs ^3] obtained from MC events. For continuum 
background, V mES ^E is the product of two orthogonal one-dimensional parametric PDFs; tties is 
well described by an empirical phase-space threshold function |14j and AE is parameterized with a 
second degree polynomial. The parameters of the continuum PDFs are floated in the fit, with rn.ES 
constrained to masses below 5.29 GeV/c 2 . ANN is described by the product of an exponential 
and a polynomial function for continuum background and by a Crystal Ball function ^5] for all 
other modes. For V m ^ n , one-dimensional non-parametric PDFs obtained from MC events are used 
to describe all modes except the signal mode itself, which is described by a non-relativistic Breit- 
Wigner line-shape. The parameters for this PDF are held fixed to the MC values and varied within 
errors to estimate systematic uncertainties. 

| Scaled p Helicity | 
160 

140 
120 
100 
80 
60 
40 
20 

°-1 -0.8 -0.6 -0.4 -0.2 -0 0.2 0.4 0.6 0.8 1 

cos 6 Hel 

Figure 1: Distribution of cos#Hei for ^-es > 5.265, \ AE\ < 0.1, and ANN > 0.85. Data are shown 
in black, with error bars. The total PDF and the B^ — » p^ir contribution are overlaid in solid 
blue and dashed red lines, respectively. 

A number of cross checks confirm that the fit is unbiased. In 1000 separate MC pseudo- 
experiments we generate the expected number of events for the various fit components before using 
the maximum likelihood fit to extract the yields and asymmetries. The distributions for each 
component are generated from the component's PDF, giving values for the fit variables mES) AE, 
ANN, and m n7T . The expected number of events is calculated from the branching fraction and 
efficiency for each individual mode. The generated number of events for each fit component is 
determined by fluctuating the expected number according to a Poisson distribution. The test is 
repeated using samples with differing asymmetry values. We repeat these MC studies using fully 
simulated signal B^ — > p^ir® events instead of generating the signal component from our PDFs. 
This verifies that the signal component is correctly modeled including correlations between the fit 
variables. As another cross check we compare the distribution of the helicity angle #Hei between 
the momenta of the charged pion and the i?-meson in the p rest frame in data with that modeled in 
MC samples for a variety of cuts. Fig. ^ shows the distribution of cos#Hei for a pseudo-signal-box 
defined by ?tt,es > 5.265, \AE\ < 0.1, and ANN > 0.8. We generally find our PDFs in good 
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agreement with the data. Finally, omitting fit variable has no significant influence on the 

signal yield, indicating that our treatment of p* background is indeed effective. 

5 Systematic Uncertainties 



Table 2: Breakdown of systematic uncertainties. 



Absolute uncertainties on yields 


Source 


a™ (Events) 


Background normalization 


+14.1 
-13.4 


PDF shapes 


+ 4.7 
- 4.2 


SCF fraction 


±12.2 


Mistag fraction 


± 2.0 


AE shift 


± 2.6 


Total 


±19 


Relative uncertainties on 


B(B ± -» ^vr ) 


Source 


4st.(%) 


Efficiency estimation 


±7.3 


B counting 


±1.1 


Total 


±7.4 


Uncertainties on Acp 


Source 


a A C p 

°Syst. 


Background normalization 


±0.006 


Background asymmetry 


±0.024 


PDF shapes 


±0.001 


Total 


±0.02 



Individual contributions to the systematic uncertainty are summarized in Table |2j 
Absolute uncertainties on yields. We calculate the uncertainty of the continuum back- 
ground estimation directly from the fit to data. The backgrounds from B decays are determined 
from simulation and fixed according to their efficiencies and branching fractions. The largest indi- 
vidual contribution comes from the B — > a^ir channel. For those individual decay modes which 
have been measured, we vary the number of events in the fit by their measured uncertainty. For all 
others we vary the amount included in the fit by ±100%. For the b — > c component, we fix the rate 
based on the number calculated from MC samples and vary the amount based on the statistical 
uncertainty of this number. The shifts in the fitted yields are calculated for each mode in turn and 
then added in quadrature to find the total systematic effect. The largest individual contribution 
comes from the p* estimation. 
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To take into account the variation of the two-dimensional non-parametric PDFs used for AE 
and m-ES) we smear the MC-generated distributions from which the PDFs are derived. This is 
effectively done by varying the kernel bandwidth ^H] up to twice its original value. For 
and ANN, the parameterizations determined from fits to MC events are varied by one standard 
deviation. The systematic uncertainties are determined using the altered PDFs and fitting to the 
final data sample. The overall shifts in the central value are taken as the size of the systematic 
uncertainty. 

We vary the SCF fraction by a conservative estimate of its relative uncertainty (±10%) and 
assign the shift in the fitted number of signal events as the systematic uncertainty of the SCF 
fraction. 

To account for differences in the neutral particle reconstruction between data and MC simu- 
lation, the signal PDF distribution in AE is offset by ±5 MeV and the data is then refitted. The 
larger of the two shifts in the central value of the yield is 2.2 events, which is taken as the systematic 
uncertainty for this effect. 

Relative uncertainties on the branching fraction. Corrections to the 7r° energy resolution 
and efficiency, determined using various control samples, add a systematic uncertainty of 7.2%. A 
relative systematic uncertainty of 1% is assumed for the pion identification. A relative systematic 
uncertainty of 0.8% on the efficiency for a single charged track is applied. Adding all the above 
contributions in quadrature gives a relative systematic uncertainty on the branching fraction of 
7.3%. Another contribution of 1.1% comes from the uncertainty on the total number of B events. 

Uncertainties on the charge asymmetry. To calculate the effects of systematic shifts in 
the charge asymmetries of background modes, each mode is varied by its measured uncertainty. 
For contributions with no measurement, we assume zero asymmetry and assign an uncertainty of 
20%, motivated by the largest charge asymmetry measured in any mode so far |TJj]. The individual 
shifts are then added in quadrature to find the total systematic uncertainty. In addition, the effect 
of altering the normalizations of the B backgrounds affects the fitted asymmetry. The size of the 
shift on the fitted Acp is taken as the size of the systematic uncertainty. 

6 Results 

The central value of the signal yield from the maximum likelihood fit is 357 ± 49 events, over 
44840 ±217 continuum events and an expected background of 872 ±62 events from other B decays. 
We find a branching fraction and charge asymmetry of 

B(B ± -> p ± 7T°) = [10.0 ± 1.4 (Stat.) ± 0.9 (Syst.)) x 10~ 6 
A C p{B ± -» p ± 7r°) = -0.01 ± 0.13 (Stat.) ± 0.02 (Syst.). 

Compared against the null hypothesis, the statistical significance ln(CNuii/^-Max) of the yield 
amounts to 8.7 standard deviations. 

The results of the fit are illustrated in Fig. |3 The plots are enhanced in signal by selecting 
only those events which exceed a threshold of 0.1 (0.05 for ANN) for the likelihood ratio R = 
(N Si ^r Sis )/(N Sis r Sig + E l ^"f kg ^f kg ), where N are the central values of the yields from the fit 
and V are the PDFs with the projected variable integrated out. This threshold is optimized by 

maximizing the ratio S = N Si % e Sis /y/N Si e e Si § + £V ivf kg ef kg where e are the efficiencies after 
the threshold is applied. The PDF components are then scaled by the appropriate e. 
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Figure 2: Likelihood projection plots for the four fit variables, (a) mES> (b) AE, (c) ANN, and (d) 
m^TT- In each plot the solid blue line represents the total PDF, the dotted green line represents the 
total background, the dotted-dashed light blue line represents the continuum contribution, and the 
dashed red line represents the signal component. The plots contain a subset of the events defined 
by a likelihood ratio of at least 0.1 (0.05 for ANN). 
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7 Conclusions 



We have measured the branching fraction and charge asymmetry for the decay — > p^n using 
a maximum likelihood fit. We obtain B{B ± -> p ± tt°) = [10.0 ± 1.4 ± 0.9] x 10~ 6 , and A C p = 
—0.01 ± 0.13 ± 0.02, respectively, where the first error is statistical and the second error systematic. 
The statistical significance of the signal is calculated to be 8.7 standard deviations. The results are 
in good agreement with the previous measurement [3]. 
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